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ABS TRACT 

Extray olation techniques including the geostrophic and gradient 
thermal wind, average wind shear from United States weather Bureau 
meximum-wind analyses, and persistence were statistically compared 
for January 1960 over the Ynited States. Calculations were made for 
layers contained between 200 mb and 10,000 or 15,000 feet, the thick- 
nesses averaging 2100 or 6500 feet, respectively. Graphic results are 
presented for both speed and direction errors for each type of extrapol- 
ation. 

it is shown that low wind speed, cross-isobar flow, and large 
curvature are associated with significant extrapolation errors. 

The author is grateful for the patience and expert counsel of 
Professor lk. J. Renard of the United States Naval rostgraduvate School, 


Monterey, California. 
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lL. ImGroduct ron 

The progress of jet aircraft to higher and higher operat- 
ing altitudes has created a demand for more precise upper- 
atmosphere information, At present, there is a particular 
need for an accurate description of the wind field in and 
about the tropopause. The meteorological parameters of most 
concern to aircraft operating at these high altitudes are the 
enroute winds and clear air turbulence. The former is normally 
obtained from soundings or other types of direct wind (or 
pressure) observations. 

Actual upper level wind reports are frequently missing 
for such areas as oceans, unpopulated continental sectors, or 
enemy held territory in wartime, This deficiency may be 
corrected by vertical extrapolation of wind from a level of 
known data by employing various aids such as the horizontal 
temperature field, empirical values of the wind shear, or 
layer of maximum wind analyses, tycume a few Although each 
of these parameters have been studied individually by others, 
no comprehensive comparison of vertical extrapolations of the 
horizontal wind at high altitudes has been published. 

This shortcoming prompted a comparative analysis of the 
following vertical extrapolation techniques in the vicinity 
of the tropopause: geostrophic thermal wind, gradient thermal 


wind, average vertical shear of the actual wind, and persistence. 





2. Previous extrapolation studies 

Various methods are presently employed to determine the 
horizontal wind field above the highest reported wind. The 
simplest extrapolation procedure considers the wind vector to 
remain constant with height. This is the best technique to 
use for short vertical distances (500 feet) or in areas where 
little change in the wind vector with height is indicated. 
Determination of the latter areas is dependent upon an 
accurate subjective analysis of the wind field below the 
desired level, on the dynamics of the wind field as associated 
with the temperature field and the location of the jet axis. 
The condition of LV faz =o is least likely to be met 
a few thousand feet above and below the tropopause and level 
of maximum wind. 

The geostrophic thermal wind is often used for extra- 
polation through layers in which the wind field is changing 
rapidly. This method requires the determination of the mean 
horizontal temperature field for the layer. Lacking a mean 
temperature field, the isotherms at some level within the layer 
are often substituted for that of the former. 

A few synoptic studies of the geostrophic thermal wind 
have been made. Ina recent survey by Cunningham [ 2 | accurate 
wind and temperature data were obtained from a specially 
instrumented B-29 flying at 280 and 383 mb over the United 
States on tracks perpendicular to the jet stream. Percent 


departures of the geostrophic and gradient thermal wind 





extrupolations from the observed shear decrease c«ponen= 

tially with increasing distance (d) over which the horizontal 
temperature gradient was measured. Results were only slightly 
better for the gradient computations. Cunningham considered the 
lack of a large difference between the two methods due to the 
small isobar curvatures in the sample (maximum curvature 
approximately 1/1000 nautical miles). 

To illustrate the advantage of using a mean wind, 
Cunningham averaged the wind over d before vectorially adding 
the geostrophic thermal wind to obtain the upper level wind. 

The departure of this wind from the actual wind varied inversely 
with d, and decreased from 11 percent at 33 miles to 2.7 

percent at 181 miles. Even better results (two percent improve- 
ment) were achieved when using a mean verification wind as 
averaged over d, rather than a spot wind. 

In order to understand the wind distribution about the 
jet core and perhaps to subjectively adjust vertical wind shear 
extrapolations, a knowledge of the jet stream is essential. 
Several empirical investigations of vertical shear in the 
vicinity of the jet stream have been recently published. 

Endlich [ 3} has developed a model cross section through the 

jet stream for use as an aid to forecasters in filling in areas 
of missing data. It provides an average wind distribution in 
a plane perpendicular to the jet in terms of percent of the 


speed at the jet core. 





In a report on the layer of maximum wind, Reiter [6] 
indicates that, among other things, there is a lack of systema- 
tic dependence of the vertical decay of the wind on the 
magnitude of the maximum wind. However, he does show percent 
vertical decay as a function of the horizontal distance from 
the jet core. A later United States Navy AROWA study [ 10) 
obtained jet stream vertical shears that showed an increase 
with jet core speed. The magnitude of vertical shears above 
the core were indicated to be as much as 30 percent greater 
than those below, However, a United States Air Force study 

{ 9| of winds over 100 knots indicates the maximum vertical 
shears to lie below the maximum wind level. 

A particularly useful description of the jet streams over 
the Pacific Ocean is provided by Serebreny (8, who presents 
numerous rules of thumb from statistical studies of these 
jets. A more comprehensive theoretical presentation of jet 


stream information is provided by Riehl et al [7] ° 





Theoretical acvclopimwat 


The equation of horizontal motion in natural coordinates can 


be written as | 5 
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and as a 


= 0, 


This becomes the gradient wind equation for V 


special case, the geostrophic wind equation when, additionally, 


In a manner similar to the development by Forsythe {4\ , 


equation (1) is differentiated with respect to p to obtain 
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equation (2) becomes 


Interchanging the Vp and 2 operators, 
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The tangential acceleration, V¥, in the first term of (7) is 
equal to the component of the pressure force in the direction 
GE the wind @nd is the mos@ydieemeult to evaluate, The first 
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term may be expanded by setting 
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COMcCOlL anasysis is questionable. Uowever;, when values of 4 
definitely greater than zero appearc@, they were taLulated as 


Ail ia@ieation of large V, 





Im this study, 2Z0@ 1m t@s been selected as the known level. 
Littic woisture exists at this height and the temperature averages 
approximately -55°C (218°K) over the United States for the period 
seleetcuetor study the temperature ata is assumed invariant 
from 200 mb to the extrapolated level with ak = Ve ne 
(temperature measured at 200 mb). With the further assumption 


that K, and V are identically zero, equation (7) becomes the 


geostrophic thermal wind equation 


OVg5 = - Vp Txk (10) 
oles £7 


expressed in its scalar working form as 
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where V;_ is obtained by applying a geostrophic wind scale for a 
contour interval of 200 feet, to isotherms at a 5°C interval. 
The gradient thermal wind equation is determined from 
equation (7) by setting V equal to zero and using the above 
Coppi CLemmonetie vector Q Vga e The equation is shown 
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Equation (12) is expressed in this form for ease in handling 
the term containing the geostrophic thermal wind vector. Figures 
lL and 2 illustrate schematic wind extrapolations using equations 
(10) and (12), respectively. 

The horizontal trajectory curvature, Kp, is obtained from 


Blaton's equation 


KK, (I+ § ov) ") 
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Figure 1. Schematic Geostrophic Thermal Wind Vector Diagram 





Figure 2, Schematic Gradient Thermal Wind Vector Diagram 





4. Computation procedures 

All data were obtained from 12002 January 1960 United States 
Weathe, Bureau facsimile charts and concurrent upper wind reports. 
20C: ab was selected as the level of known Gata because, averaging 
39,000 feet for this data sample, it was considered to be suf- 
ficiently close to both the level of maximum wind and the present 
jet aircraft maximum operating altitude. At each station extra- 
polations of winds to 40,000 and 45,000 feet were compared to 
reported winds at those levels. United States Weather Bureau 
tropopause-height and maximum-wind analyses oe this period were 
also utitizee. | 

a. Thermal Wind extrapolations 

Stations were chosen where a well defined thermal zvradient 
existed from which the wind factor, Ves could be measured with a 
Peostropiie wing Scale acot@eqcer teu si: Bection.3. Vy_ Ls multi-= 

i ac 

plied by .114z2 per thousand feet pomevtain (A) in equation (12) 
and Figure l. V, is the vector wind at 200 mb . The hone at 


the upper level (Vz) is obtained by the vector (graphical) 


addition of A Vs to V1 > 


For term (B), isobar curvature (K,) is measured by the use 
of overlay arcs corrected for map scale variation. The 200-mb 
wind (V,) and associated pressure-system movement (c) are 
determined and entered in equation (13) to yield trajectory 
curvature, K,, considered constant through the layer. 


The scalar Vgr is approximately equal te V, the actual 


Meanewiid Met the layer. As a first estimate of V , an 


10 





approximate vector mean ot Vg and Vig may be obtaineu. An improve- 
ment is provided by the#vector mean of Vr and V1: 

Term (A), /Vus, is useu for the first estimate of AVar. 
The coriolis parameter (f£) is taken from published tables (i | ‘ 
Although (B) is a vector parallel to INGE. it may be drawn 
parallel to /\Vzs as a first approximation. 

In term (C) AVgr is initially considered equal to the 
scalar magnitude of /\Vgs. Vector (C) parallels Ver, approxi- 
mately by un which may be obtained as previously described or by 
trial and error. 

The resultant Ver is reentered in (B) and (C), Var is 
recomputed and a second approximation to A\Vgr calculated. This 
is generally sufficient to produce the desired accuracy in 

2 Gee 

b. ShearewWini extrapolations 

Over the United States, the facsimile maximum-wind anal- 
yses provide the height and velocity of the maximum wind and 
the vertical wind shear in the layer containing the level of 
maximum wind. (This level is hereafter referred to as the 
LMW). These shears are a subjectively determined average of 
the plotted shears 10,000 feet above and below the LMW. 

The shear obtained from the arene is multiplied by the 
thickness of the extrapolation layer (in thousands of feet). This 
gives [\Vs which is added to or subtrected from the wind at the 
known level (200 mb), depending upon whether it is below or 


above the LMW. For extrapolations crossing the LMW, a net /\Vs 


ll 





is obtained from the difference between . Vs from 200 mb to 
the LMW and —.Vs from the LMW to the desired upper level. The 


wind direction is considered to remain constant through the 


layer. 


A shearewind computation was made for all stations for 


which a thermal wind computation was made, 


i 





5. Graphic analysis 901 Chowputations 

The objective of this study was to statistically consider 
the error distribution of various methods of vertical extra- 
polation of the horizontal wind field and its relation to 
other parameters, and thus determine the most accurate extra- 
polation procedure. A graphic representation was selected as 
best meeting these requirements. 

a. Error discussion 
Error as used herein is defined as: 
(1) the angle between the extrapolated wind direction and 
the actual wind direction (from teletype reports) at a given 
level. A plus error means that the extrapolated wind had 
turned through a greater angle (and in the same direction) than 
the actual wind for the given layer, the bottom of which is at 
200 mb. A minus error represents too little turning (or turning 
in the wrong direction) of the extrapolated wind. Therefore, 
for extrapolations uSing persistence, all errors are minus. 
However, since the distributions of the thermal. wind~direction 
extrapolations are nearly normal, a clearer comparison is 
obtained by distributing the persistence extrapolation errors 
symmetrically about an error of zero, 1. e. for any given error, 
half the cases are considered plus, half are taken as minus. 
(2) the difference between the extrapolated wind speed and the 
actual wind speed (from teletype reports) divided by the actual 
wind speed at the level. Thus, the result is expressed in 


percent of the actual Wind speed. 
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It wWust be ceatisazel] that, at times, extreme errors may 
arise solely due to errors in observing and reporting the winds. 
Reiter _ 6) has shown how the method of determining upper winds 
can pruduce the unrealistic oscillations sometimes found in 
vertical wind profiles. He also described how the abbreviated 
teletype codes can remarkably alter the actual upper~level 
sounding curve. His illustrations particularly emphasize the 
desirability of averaging and smoothing to obtain representative 
wind profiles. 

Another error arises when the neglected term, the 
tanzential acceleration becomes large. Subject to errors in 
analysis, cross-isobar flow is an indication that this accelera- 
tion is significant. To determine the effect of this term, 
the number of clearly apparent cross-isobar angles greater than 
15 degrees have been recorded in boxes (()) on the graphs which 
follow. 

The number of winds under 20 faate have also been recorded 
in circles ((> ) on the graphs. At this low wind speed, 
large errors in direction and percent speed result from small 
wectOr Grrors:. 

b. Data 

Information from 150 soundings were utilized to obtain 
107 layers whose thickness averaged 2100 feet (range from 
900 to 3300 feet) and 149 layers whose thickness averaged 6500 
feet (range from 3300 to 8400 feet). Two thermal wind extra~ 
polations and a shear extrapolation were compared with persis-~ 


tence for each layer, considering direction and speed separately. 
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Of the 150 reperts, 21 showed a near-zero therval pradient 


at the 200-mb level, 45 bad trajectories which were nearly 


Straight (curvature less than 1/2500 nautical miles), 22 cases 


were associated with anticyclonic curvature, and 62 with 


cyclonic curvature. 


The following factors were considered for their effect on 


the extrapolation crror: 


|e 


2. 


Cross-isobar angle 


2 


Location of the jet maximum in relation to the station 
(horizontal distance and direction) 


Tropopause height 
Level of maximum wind 


Whether wind speed was increasing or decreasing with 
height 


Angle between the thermal wind and actual wind 
Magnitude of wind speed 
Magnitude and algebraic sign of trajectory curvature 


Magnitude of the thermal gradient 


The graphs that follow represent the most significant results 


of the above considerations. 


Figure 3. 


This graph presents frequency of errors in the extra- 


polated angular turning of the wind for layers averaging 2100 


feet. Little deviation from persistence is shown by the geo- 


Strophic or gradient thermal wind errors, since 87 percent of 


the geostrophic or gradient thermal wind computations lie withiu 


ae : 
- 15 degrees of the actual wind. Of the persistence cases, 


JES 





84 percent tall within |f\ie range. 

The slight skew toward negative an;,le error wiaitates that 
the cwaputations, on the average, uncerestimate the actual 
angular turning. A further analysis of these winds indicates 
that this slight skewness is due to cases in which the wind 
speed decreases with increasing height. 

Considering only the cases of Kj,>1/2000 nautical miles, 
where Q Vgs exceeds two knots per thousand feet, does not 
change the substance of the foregoing conclusions. 

As would be expected for extreme angular errors 
(> | 40 degrees| ), there are relatively large percentages of 


low-wind-speed, cross-isobar and large-curvature cases. 


Figure 4 

This figure presents the ears type of information displayed 
in Figure 3, but only for extrapolations through layers averig- 
ing 6500 feet. For =e degrees and as degrees errors, 
respectively, the included percentage frequencies are: 63 and 
76 percent for persistence, 70 and 84 percent for geostrophic 
and 73 and 88 percent for gradient thermal wind extrapola- 
tions. The latter represents a significant improvement over 


persistence, although slighly skewed. 
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The sevcer@ect ske"iog is cbserved for the geostrophic cases, 
iv, ceumtrast to figure 3. However, again, this effect ia dug to 
the cases of winds which show decrease with increasing height. 

At the extremes of the distribution, low wind speeds account for 
a large portion of the error. 

Large curvature is present in most of the cases associated 
with errors greater than 20 degrees. Note that the mode for 
cases of large curvature lies at -10 degrees. Even so, 74 
percent of the geostrophic> and 82 percent of the gradient-thermal 


: : : : : AA 
wind cases in this category lie within - 25 degrees error. 
Figure 5 
The distribution of the percent error in the extrapolated 
wind speed from the actual (reported) wind speed is shown here, 
for layers averaging 2100 feet. All four types of computations 
| | ; ; , AM. 
peak together at zero with Little skew cxeept possibly the AY¥s 


ir 


Pees 2 . 
wind cases. Within - 15 percent and - 25 percent, respectively 


lics &6 and 60 percent of the persistence cases, 66 and 85 per~ 
cent of the shear cases, /1 and 89 percent of the geostrophic - 
and 75 and 90 percent of the gradient-thermal wind cases. 
Persistence and snear show the maximum large error samples. 
As expected, they are associated with a correspondingly large 
number of cross-isobar and low-wind-speed cases, Cases with 
large curvature are considered on the next two graphs. 
Figure c 
This figure depicts the error distributions for those 


extrapolations in figure 6 which are associated with winds 
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decreasing speed with increasing height. The curves in figures 
5 gnd 6 are similar, cxcept that persistence occurs only on 

the positive erroz side (because winds are decreasing 

Figure o shows somewhat more irregularity than figure 5 because 
of the smaller sample in the former. Presuming the plus skew 
of the geostrophic thermal wind extrapolation curve is valid, 
it indicates the non-applicability of the simple geostrophic 
thermal wind for large curvature cases ~ a not unexpected 
result from theory in view of the high percentages of cyclonic 
cases consicered here. 

Vigure 7 

Fizure 7 depicts the frequency of speed errors for 
eases of winds increasing with increasing height. Here is 
shown a pronounced skew, with a mode at -10 percent, for all 
techniques except persistence. The number of large curvature 
cases is again too small to be of great significance, although 
the curves follow those of figure 5 quite closely. 

The overall skew indicates that none of the extrapola- 
tion techniques increases the wind speed adequate ly through 
the iayer, the computations being on the average ten 
percent too low. 

Figure 8 

This picture indicates the extreme variability in the 
accuracy of wind speeds extrapolated through a layer averag- 
ing 6500 fect. Thermal wind computations show major peaks 


at both zero and -20 percent error so that a further error 
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investigation was made necessary (figures 9 and 10). Persistence 


shows a broader curve, peaking at +10 percent and is apparently 
of less significance. Shear computations show a rather irregular 
distribution, with the mean approaching -20 percent error. 

Large curvature cases have a mode at +10 percent with a 


secondary peak at -20 percent error. 
Figure 9 

Considering only the cases of wind which decreases speed 
with increasing height from figure 8, a much clearer mode at 
+10 percent error is obtained. These curves indigate that the 
computations either underestimate the decrease of wind with 
height or increase the wind on the aveyage, by about ten 
percent. 

The data for the secondary maximum at -20 percent error 
was examined, and in addition to several low-wind-speed and 
cross-isobar cases, there are four that cross the tropopause 
and six that pass through the layer of maximum wind. Some of 
the explanations for the occurrence of this type of error would 
be a secondary maximum above the primary maximum in the actual 
wind profile, an excessive thermal gradient at the level where 
the extrapolation began, or a thermal gradient associated with 
a relative minimum just above the original level. 

The error curves indicate that the best results could be 
obtained by using persistence, if the wind was known to be 
decreasing. <A wind at the top of the layer correct within 
+ 20 percent would be obtained 60 percent of the time by sub- 
tracting 15 percent from the wind at the known level. The 
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extrdfe ervors {> 00 percent) that may arise, are iu part attri- 
butaLle to tue laryc number of cross-isobar and low-wind-speed 
cases Unere. 

The curve for large curvature cases follows the parent 
curve in figure 8 except for the indication that here the geo- 
strophic thermal wind extrapolations do not decrease the wind 
adequately. 

Pooune. 107 

For the cases of winds increasing with increasing height 
in figure 8, shear extrapolations produce a maximum skew with 
a mode of ~30 percent error. Both thermal wind computations 
have modes of zero, with 60 percent of the geostrophic and 66 
percent of the gradient thermal wind computations lying in the 
range -5 20 percent. 

Persistence would be better, however, when the wind 
speed is known to be increasing, as 78 percent of the winds 
lie within 20 percent of a meal “set =lompereene. 

Large-curvature cases do not show any clear indications of 
relevant trends, 

Hie over 

The results in figures 5 - 10 indicate that it is 
important to know whether the winds are increasing or dcecreas- 
ing through a layer. Table 1 was assembled to predict the 
accuracy with which each extrapolation technique could give 
this iniormation. Except for the shear computation through 


the thicker layer, the correct result can be expected 60 per- 


cent of the time. 
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Cea conclusions 

From-this study, it iS@apparent= that more accuraterextra- 
polations in the vicinity of 200 mb can be made if the influence 
of various restrictive condit#ins is also considered. A 
qualitative estimate of the likelihood of encountering ieee 
errors is provided by the proportionately large number of 
cross-isobar and low-wind-speed (and possibly, large~curvature) 
cases occurring near the limits of the error distribution 
graphs. In addition, the type of extrapolation produces an 
intrinsic error. Even more important is the determination 
of whether the wind speed is increasing or decreasing 
with height, 

Considering wind direction alone, for a 2100-foot layer, 
PersistGmeemus observed Lo be aceveod aseciimemeen, Ehe thermal 
wind extrapolations. For a 6500-foot layer the gradient -" 
thermal wind computation is superior, with 73 percent falling ¥ 
within = 15 degrees and 88 percent within : Zomdesrees of 
actual turning through the layer. For large curvature cases 
in a 6500-foot layer, both thermal wind computations produce 
a small error that averages five degrees too little turning. 
wherd computing time is limited, pe renocenees for these same 
cases, will produce an error less than 25 degrees 76 percent 
of the time. 

Better results are obtained for vertical extrapolations 
of the speed using thermal wind extrapolations (thermal 


se + 
e@rauient accurate to - 15 percent and ~- 25 percent 75 and 
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JO percent of the cipa ecopectively) compareu to persistence or 
sheat Gowputations Ler a 2100-foot layer. For uecreasing wine 
spceu Vith height, geostrophie thermal wind computations do not 
decreasc the wind sufficiently, whereas for increasing speed 
with height, all computations average ten percent low. 

Extrapolation results for speed are least accurate for a 
6500-toot layer. Generally, best results are obtained for 
shear computations with the speed error in the range -10 ‘ Jaa) 
percent, 70 percent of the time. For decreasing wind speed 
with increasing height, all curves peak near plus ten percent 
error with a secondary mode near -20 percent. Consequently, 
persistence minus 15 percent. is correct within 226 percent 
of the true wind, 60 percent of the time. Even sa, errors 
greater than 50 percent occur 20 percent of the time. Consider- 
ing only increasing winds with height, thermal wind computa- 
tions peak at zero. Again, persistence plus 15 percent yields 
78 percent of the cases within E20 percent. 

ALL thermal wind extrapolations predict the correct 
direction of shear(plus or minus) in the vicinity of 200 mb 
60 percent of time as compared to 50 percent for shear.wind 
computations. 

Other approaches that are recommended for further study, 
to obtain more accurate wind extrapolations, are listed below. 
(1) Instead of considering speed and direction separately, 
consider the magnitude (and perhaps direction) of the vector 


error (especially for low wind speeds). 
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3) Ein@ cross-isobar cases separatcly. 

ey) @t’rmkae a more accurate method of predicting whether 


Wé Wina is increasing or decreasing with increasing height. 
(:) oscperatcly determine the shear above and below the LMI 
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